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Q\ . Abstract. 

In this paper we study symmetry reductions and exact solutions of the shallow water wave 
(SWW) equation 

Uxxxt + au x u xt + I3u t u xx - u xt - u xx = 0, (1) 

where a and (5 are arbitrary, nonzero, constants, which is derivable using the so-called Boussinesq 
approximation. Two special cases of this equation, or the equivalent nonlocal equation obtained 
by setting u x = U, have been discussed in the literature. The case a = 2/3 was discussed by 
Ablowitz, Kaup, Newell and Segur [Stud. Appl. Math., 53 (1974) 249], who showed that this case 
was solvable by inverse scattering through a second order linear problem. This case and the case 
a = (3 were studied by Hirota and Satsuma [J. Phys. Soc. Japan, 40 (1976) 611] using Hirota's 
bi-linear technique. Further the case a = (3 is solvable by inverse scattering through a third order 
linear problem. 

In this paper a catalogue of symmetry reductions is obtained using the classical Lie method and 
the nonclassical method due to Bluman and Cole [J. Math. Mech., 18 (1969) 1025]. The classical 
Lie method yields symmetry reductions of (1) expressible in terms of the first, third and fifth Pain- 
leve transcendents and Weierstrass elliptic functions. The nonclassical method yields a plethora 
of exact solutions of (1) with a = j3 which possess a rich variety of qualitative behaviours. These 
solutions all like a two-soliton solution for t < but differ radically for t > and may be viewed as a 
nonlinear superposition of two solitons, one travelling to the left with arbitrary speed and the other 
to the right with equal and opposite speed. These families of solutions have important implications 
with regard to the numerical analysis of SWW and suggests that solving (1) numerically could pose 
some fundamental difficulties. In particular, one would not be able to distinguish the solutions in 
an initial value problem since an exponentially small change in the initial conditions can result in 
completely different qualitative behaviours. 

We compare the two-soliton solutions obtained using the nonclassical method to those obtained 
using the singular manifold method and Hirota's bi-linear method. 

Further, we show that there is an analogous nonlinear superposition of solutions for two 2 + 1- 
dimensional generalisations of the SWW equation (1) with a = (3. This yields solutions expressible 
as the sum of two solutions of the Korteweg-de Vries equation. 
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1 Introduction. 

In this paper we discuss symmetry reductions and exact solutions for the shallow water wave 
(SWW) equation 

A = u xxxt + au x u xt + I3u t u xx - u xt - u xx = 0, (1.1) 

where a and (3 are arbitrary, nonzero, constants, which can be derived from the classical shallow 
water theory in the Boussinesq approximation (cf., [26]). This equation may be written in the 
nonlocal form (set u x = U) 

U xxt + aUU t - (3U x d x x XJ t -U t -U x =0, (1.2) 

where (3" 1 /) (x) = f°° f(y) dy; in fact this is the form of the physical model derived in [26]. 

Two special cases of (1.2) have attracted some attention in the literature, namely a = 2(3 and 
a = (3. In their seminal paper on soliton theory, Ablowitz, Kaup, Newell and Segur [2] show that 

U xxt + 2(3VV t - &U x d~ x V t -U t -U x = 0, (1.3) 

which is (1.2) with a = 2(3 is solvable by inverse scattering. Further, Ablowitz et al [2] remark that 
(1.3) reduces in the long wave, small amplitude limit to the celebrated Korteweg-de Vries (KdV) 
equation 

u t + u xxx + 6uu x = 0, (1.4) 

which was the first equation to be solved by inverse scattering [27], and they also comment that 
(1.3) also has the desirable properties of the regularized long wave (RLW) equation [8,47] 

Uxxt + uu x - u t - u x = 0, (1.5) 

sometimes known as the Benjamin-Bona-Mahoney equation, in that it responds feebly to short 
waves. Additionally, we note that (1.5) and (1.3) have the same linear dispersion relation 
u>(k) = —k/(l + k 2 ) for the complex exponential u(x,t) ~ exp{i[kx + Lo(k)t]} . However, in contrast 
to (1.3), the RLW equation (1.5) is thought not to be solvable by inverse scattering (cf., [42]). 
Hirota and Satsuma [34] studied both (1.3) and 

U xxt + pUU t - (3Uxd x l V t -U t -U x = 0, (1.6) 

i.e. (1.2) with a = 2/3, using Hirota's bi-linear method [32] and obtained iV-soliton solutions for 
both equations (see also [43]). 

In the sequel, we shall refer to (1.1) with a = 2(5 

u xxxt + 2(3u x u xt + (3u t u xx - u xt - u xx = 0, (1.7) 

as the SWW-AKNS equation and (1.1) with a = (3 

u xxxt + (3u x u xt + (3u t u xx - u xt - u xx = 0, (1.8) 

as the SWW-HS equation. 

The scattering problem for the SWW-AKNS equation (1.7) is the second order problem [2] 

^xx + \(iu x i) = Xtp, (1.9) 

with associated time-dependence 

(4A - 1)Y> 4 = (1 - (3u t )ip x + \(3u xt i>, (1.10) 
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where A is the constant eigenvalue, and ip X xt = iptxx if and only if u satisfies (1.7) We note that 
(1.9) is the time-independent Schrodinger equation which is also the scattering problem for the 
KdV equation (1.4) [27]. In contrast, the scattering problem for the SWW-HS equation (1.8) is the 
third order problem [23] 

Ipxxx + {\@u x -l)ip x = Xtp, (1.11) 

with associated time-dependence 

3AV>t = (1 - Pu t )ip xx + f3u xt i/j x . (1.12) 
We remark that (1.11) is similar to the scattering problem 

ipxxx + K 1 + 6u)V>x + \[u x - iV3d-\v t )\ V = (1.13) 
which is the scattering problem for the Boussinesq equation 

u xxxx + 3(u 2 ) xx + u xx = u t t, (1-14) 

and which has been comprehensively studied by Deift, Tomei and Trubowitz [24]. 

The SWW equation (1.1) was also discussed by Hietarinta [31] who shows that it can be 
expressed in Hirota's bi-linear form [32] if and only if either (i), a = f3, when it reduces to (1.8), 
or (ii), a = 2/3, when it reduces to (1.7). Further, in [19] it is shown that the SWW equation (1.1) 
satisfies the necessary conditions of the Painleve tests due to Ablowitz, Ramani and Segur [3,4] and 
Weiss, Tabor and Carnevale [55] to be completely integrable if and only if either a = (3 or a = 2(3. 
These results strongly suggest that the SWW equation (1.1) is completely integrable if and only if 
it has one of the two special forms (1.7) or (1.8), which are both known to be solvable by inverse 
scattering. 

The classical method for finding symmetry reductions of partial differential equations (PDEs) is 
the Lie group method of infinitesimal transformations (cf., [10,44]). Though this method is entirely 
algorithmic, it often involves a large amount of tedious algebra and auxiliary calculations which 
can become virtually unmanageable if attempted manually, and so symbolic manipulation programs 
have been developed, for example in macsyma, maple, mathematica, mumath and reduce, to 
facilitate the calculations. An excellent survey of the different packages presently available and a 
discussion of their strengths and applications is given by Hereman [30] . 

In recent years the nonclassical method due Bluman and Cole [9] (in the sequel referred to 
as the nonclassical method) and the direct method of Clarkson and Kruskal [17] have been used 
to generate many new symmetry reductions and exact solutions for several physically significant 
PDEs that are not obtainable using the classical Lie method, which represents important progress 
(cf., [16,29] and references therein). 

Symmetry groups and associated reductions and exact solutions have several different important 
applications in the context of differential equations (see, for example, [10,16,44] for further details 
and references): 

• Derive new solutions from old solutions . Applying the symmetry group of a differential 
equation to a known solution yields a family of new solutions. Quite often interesting solutions 
can be obtained from trivial ones. 

• Integration of ODEs . Symmetry groups of ODEs can be used to reduce the order of the 
equation, for example to reduce a second order equation to first order. 

• Reductions of PDEs . Symmetry groups of PDEs are used to reduce the total number of 
dependent and independent variables, for example from a PDE with two independent and 
one dependent variables to an ODE. 

• Linearisation of PDEs . Symmetries groups can be used to discover whether or not a PDE can 
be linearised and to construct an explicit linearisation when one exists. 
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• Classification of equations . Symmetry groups can be used to classify differential equations into 
equivalence classes and to choose simple representatives of such classes. 

• Asymptotics of solutions of PDEs . It is known that as solutions of PDEs asymptotically tend to 
solutions of lower-dimensional equations obtained by symmetry reduction, some of these special 
solutions will illustrate important physical phenomena. In particular, exact solutions arising 
from symmetry methods can often be effectively used to study properties such as asymptotics 
and "blow-up". 

• Numerical methods and testing computer coding . Symmetry groups and exact solutions of 
physically relevant PDEs are used in the design, testing and evaluation of numerical algorithms; 
these solutions provide an important practical check on the accuracy and reliability of such 
integrators. 

• Conservation Laws . The application of symmetries to conservation laws dates back to the 
work of Noether who proved the remarkable result that for systems arising from a variational 
principle, every conservation law of the system comes from a corresponding symmetry property. 

• Further Applications . There are several other important applications of symmetry groups 
including bifurcation theory, control theory, special function theory, boundary value problems 
and free boundary problems. 

The method used to find solutions of the determining equations for the infinitesimals in both the 
classical and nonclassical case is that of Differential Grobner Bases (DGBs), defined to be a basis B 
of the differential ideal generated by the system such that every member of the ideal pseudo-reduces 
to zero with respect to B. This method provides a systematic framework for finding integr ability 
and compatibility conditions of an overdetermined system of PDEs. It avoids the problems of 
infinite loops in reduction processes, and yields, as far as is currently possible, a "triangulation" of 
the system from which the solution set can be derived more easily (cf., [18,41,48,49]). In a sense, 
a DGB provides the maximum amount of information possible using elementary differential and 
algebraic processes in a finite time. 

The triangulations of the systems of determining equations for infinitesimals arising in the 
classical and nonclassical methods in this article were all performed using the maple package 
diffgrob2 [40]. This package was written specifically to handle fully nonlinear equations of 
polynomial type; packages such as those in [50,51,53] have been developed for the study of linear 
equations. All calculations are strictly "polynomial", that is, there is no division. Implemented 
there are the Kolchin-Ritt algorithm, the differential analogue of Buchberger's algorithm [14] using 
pseudo-reduction instead of reduction, and extra algorithms needed to calculate a DGB (as far as 
possible using the current theory), for those cases where the Kolchin-Ritt algorithm is not sufficient 
[41]. Designed to be used interactively as well as algorithmically, the package has proved useful for 
solving some fully nonlinear systems. As yet, however, algorithmic methods for finding the most 
efficient orderings, the best method of choosing the sequence of pairs to be cross-differentiated, for 
deciding when to integrate and read off coefficients of independent functions in one of the variables, 
for finding the best change of coordinates, and so on, are still the subject of much investigation. 

In §2 we find the classical Lie group of symmetries and associated reductions of (1.1), which 
are expressible in terms of the first, third and fifth Painleve transcendents and Weierstrass elliptic 
functions. Then in §3 we discuss the nonclassical symmetries and reductions of (1.1). In particular, 
the nonclassical symmetry reductions obtained for (1.8) generate a wide variety of interesting exact 
analytical solutions of the equations which we plot using maple. In §4 we compare the two-soliton 
solutions generated in §3 with those obtained using the singular manifold method and Hirota's 
bi-linear method. In §5 we show that there are analogous symmetry reductions and exact solutions 
also occurs for two 2 + 1-dimensional generalisations of the SWW equation (1) with a = (3. We 
discuss our results in §6. 
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2 Classical Symmetry Reductions of the SWW Equation. 

To apply the classical method to the SWW equation (1.1) we consider the one-parameter Lie 
group of infinitesimal transformations in (x, t, u) given by 

x = x + e£(x,t,u) + 0(e 2 ), (2.1a) 
t = t + er{x,t,u) +0{e 2 ), (2.1b) 
u = u + e<f>(x,t,u) +0(e 2 ), (2.1c) 

where e is the group parameter. Then one requires that this transformation leaves invariant the 
set 

S A = {u(x,t) : A = 0}, (2.2) 

of solutions of (1.1). This yields an overdetermined, linear system of equations for the infinitesimals 
£(x,t, u), r(x,t,u) and <f)(x,t,u). The associated Lie algebra of infinitesimal symmetries is the set 
of vector fields of the form 

v = t, u)d x + t(x, t, u)d t + <j){x, t, u)d u . (2.3) 

Having determined the infinitesimals, the symmetry variables are found by solving the characteristic 
equations 

dx = dt = du 
£(x,t,u) r(x,t,u) (f)(x,t,u) , 

which is equivalent to solving the invariant surface condition 

ip = £(x, t, u)u x + t(x, t, u)u t — 4>(x, t, u) = 0. (2.5) 
The set 5a is invariant under the transformation (2.1) provided that 

pr( 4 )v(A) =0, (2.6) 

A=0 

where pr^ 4 ) v is the fourth prolongation of the vector field (2.3), which is given explicitly in terms 
of £, r and d> (cf., [44]). This yields a system of fourteen determining equations, as calculated 
using the MACSYMA package symmgrp.max [15]. A triangulation or standard form [18,48] of these 
determining equations is the following system of eight equations, 

£u =0, £t = 0, £ xx = 0, t u = 0, t x = 0, 
a4> x - 2£ x = 0, P4> t - r t - 4> x = 0, 4> u + £ x = 0, 

from which we easily obtain the following infinitesimals, 

/ 2x t\ git) . . 

£ = kix + k 2 , T = g(t), (p = -K 1 [u +— 5- + K3, (2.7) 

where g(t) is an arbitrary differentiable function and «i, K2 and K3 are arbitrary constants. The 
associated Lie algebra is spanned by the vector fields 

vi = xd x - (u- ^ - -Q d u , v 2 = d x , v 3 = d u , v 4 (g) = g(t) (d t + f3~ l d u ) . 

Solving (2.4) with £, r and 4> given by (2.7), or equivalently solving (2.5), we obtain the following 
two canonical symmetry reductions. 
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Case 2.1 Ki ^ 0. In this case we set g(t) = f(t)/f(t), Ki = 1 and n 2 = K3 = (with / = el//dt). 
Hence we obtain the symmetry reduction 

u(x, t) = f(t)w(z) + - + i z = xf(t), (2.8) 
a p 

where w(z) satisfies 

zw"" + Aw'" + (q + 0)zw'w" + f3ww" + 2q (u/) 2 = 0, (2.9) 

where ' = d/dz. It is straightforward to show using the algorithm of Ablowitz et al [4] that this 
equation is of Painleve-type, i.e., its solutions have no movable singularities other than poles, only 
if either (i), a = (3 or (ii), a = 2(3. These two special cases (2.9) are solvable in terms of solutions 
of the third Painleve equation [36] 

d 2 y 1 / dy\ 2 1 dy „ by 2 + c d 

' ' ' y +ay 3 + — + -, (2.10) 



dx 2 y \dx J x dx x y 

and the fifth Painleve equation, 

d 2 y ( 1 1 \ / dy V 1 dy (y - l) 2 / 6\ cy dy{y + 1 



dx 2 \2y y — 1 ) \dx ) x dx x 2 \^~^y)~^x^~ y — 1 ' ^ ^ 

with a, 6, c and d constants (see [19] for details). Hence the Painleve Conjecture [3,4] predicts that 
a necessary condition for (1.1) to be completely integrable is that (a — (3)(a — 2(3) = 0, i.e., only if 
(1.1) has one of the two special forms (1.8) or (1.7). 

Case 2.2 n\ = 0. In this case we set g(t) = l/f(t), k 2 = 1 and K3 = —1/(3 and obtain the 
symmetry reduction 

u(x,t) =w(z) +t/(3, z = x-f(t), (2.12) 

where W(z) = w'(z) satisfies 

(W'f + ±(a + (3)W 3 = AW + B, (2.13) 

with A and B arbitrary constants. This equation is equivalent to the Weierstrass elliptic function 
equation 

(p'f =4p 3 -g 2 p-g 3 , (2.14) 

where g 2 and 53 are arbitrary constants [56] . 

We remark that the vector field v 4 (g) shows that (1.1) is invariant under the following variable 
coefficient transformation 

x = x, t = g(t), u = u+[g(t)-t}/(3, (2.15) 

i.e., if u(x, t) is a solution of (1.1), then so is u(x, t). Further, the associated Lie algebra q = {v 4 (<7)}, 
with g(t) G C°°(E), is a Virasoro algebra since 

[v 4 (yi),v 4 (y 2 )] = v 4 {g 1 g 2 - g 2 g x ) , 

where [A, B] = AB — BA is the standard commutator. Virasoro algebras frequently arise in the 
study of symmetry reductions for 2 + 1-dimensional equations such as the Kadomtsev-Petviashvili 
and Davey-Stewartson equations (see, for example, [57] and the references therein), but not 1 + 1- 
dimensional equations such as (1.1). 
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It might be thought that the presence of the arbitrary function g(t) in the infinitesimals 
associated with the SWW equation (1.1) is an artefact of the fact that we are considering (1.1) 
rather than (1.2), which is the physical equation. In view of the nonlocal term, we write the 
nonlocal equation (1.2) as the system 

U xxt + aUU t + 0VU x - U t - U x = 0, V X = U U (2.16) 

where a and (5 are arbitrary, nonzero constants. To apply the classical method to (2.16), we consider 
the one-parameter Lie group of infinitesimal transformations in (x, y, U, V) given by 

i = x + e^(x,y,U,V) + 0( y e 2 ), 
y = y + e^ 2 (x,y,U,V) + 0(e 2 ), 
U = U + eMx,y,U,V) + 0(e 2 ), 
V = V + eMx,y,U,V) + 0(e 2 ), 

where e is the group parameter. Solving the resulting determining equations yields the infinitesimals 

£ = kix + k 2 , r = g(t), 0! = -2m(U- 1/a), 2 = - (ki + ^ (V - 1/(3), 

where K\ and k 2 are arbitrary constants and g(t) is an arbitrary differentiable function and the 
associated vector fields are 

w 1 = xd x -2(U-l/a)d u + (v-l/(3)d v , w 2 = d x , w 3 (g) = g(t)d t + ^ (V - \/(3)d v . 

This means that the system (2.16) is invariant under the transformation 

Hence if {U(x,t),V(x,t)} are solutions of (2.16), then so are {U(x,i), V(x,i)}. 

3 Nonclassical Symmetry Reductions of the SWW Equation. 

There have been several generalisations of the classical Lie group method for symmetry 
reductions. Bluman and Cole [9], in their study of symmetry reductions of the linear heat equation, 
proposed the so-called nonclassical method of group-invariant solutions. This method involves 
considerably more algebra and associated calculations than the classical Lie method. In fact, it 
has been suggested that for some PDEs, the calculation of these nonclassical reductions might be 
too difficult to do explicitly [45], especially if attempted manually since the associated determining 
equations are now an overdetermined, nonlinear system. For some equations such as the KdV 
equation (1.4), the nonclassical method does not yield any additional symmetry reductions to 
those obtained using the classical Lie method, while there are PDEs which do possess symmetry 
reductions not obtainable using the classical Lie group method. It should be emphasised that the 
associated vector fields arising from the nonclassical method do not form a vector space, still less 
a Lie algebra, since the invariant surface condition (2.5) depends upon the particular reduction. 
Subsequently, these methods were further generalised by Olver and Rosenau [45,46] to include 
"weak symmetries" and, even more generally, "side conditions" or "differential constraints", and 
they concluded that "the unifying theme behind finding special solutions of PDEs is not, as is 
commonly supposed, group theory, but rather the more analytic subject of overdetermined systems 
of PDEs" . 
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In the nonclassical method one requires only the subset of Sa given by 

Sa,4, = MM) = A(u) = 0, VH = 0} , 



(3.1) 



where S/\ is as defined in (2.2) and ^ = is the invariant surface condition (2.5), is invariant under 
the transformation (2.1). The usual method of applying the nonclassical method (e.g., as described 
in [39]), to the SWW equation (1.1) involves applying the prolongation pr^ 4 ) v to the system of 
equations given by (1.1) and the invariant surface condition (2.5) and requiring that the resulting 
expressions vanish for u G <Sa,v>j i- e -> 



pr< 4 > v (A) 



A=0,i/>=0 



0, pr (1) v(V>) 



0. 



A=0,V>=0 



(3.2) 



It is easily shown that 



pr (1) v(ip) = - {i u u x + T u u t - <pu) ip, 



which vanishes identically when ip = without imposing any conditions upon £, r and <fi. However 
as shown in [20], this procedure for applying the nonclassical method can create difficulties, 
particularly when implemented in symbolic manipulation programs. These difficulties often arise 
for equations such as (1.1) which require the use of differential consequences of the invariant 
surface condition (2.5). In [20] we proposed an algorithm for calculating the determining 
equations associated with the nonclassical method which avoids many of the difficulties commonly 
encountered; we use that algorithm here. 

In the canonical case when r ^ we set r = 1. (We shall omit the special case r = 0; in that 
case one obtains a single condition for 4> with 424 summands, and even considering the subcase 
<j) u = leads to an equation more complex than the one we are studying.) Eliminating ut, u x t and 
^xxxti HI (1.1) using the invariant surface condition (2.5) yields 

A = 4 ) uV'xxx ~i~ 3<f> uu U x U xx + 3<f> ux U xx + 4 ) uuu'U' X ~\~ 3(f) xuu U x + 3(f) xxu U x + 



%xxx^x 3^ xx 1l^ 
,2 



3^ X U XXX 



3& 



,.u" 



xxu^x i Yxxx 
x^x.r 



Q£,uuU x U x 



3£, XX u^ J x ^Cxu^x^xx £,u i^^x^xxx ~t~ ^xi) C^xxxx 

+ (au x - 1) [4> x + <j) u u x - £ x u x - £ u u\ - £u xx \ + u xx [(3 {4> - £,u x ) - 1] = 0, (3.3) 

with t a parameter, which involves the infinitesimals £ and 4> that are to be determined. Now we 
apply the classical Lie algorithm to this equation using the fourth prolongation pr^ 4 - 1 v and eliminate 
^^xxxx USUI g (3.3). This yields the following overdetermined, nonlinear system of equations for £ 
and <fi (contrast the classical case discussed in §2 above where they are linear). 



4>uu = o, 

(a + /?)(&,+ Cx) = 0, 

£4>tu + 3£ 2 £ xx + 3£ t £r - 3££ xt - 3£, 2 4>xu - £t4>u = 0, 

a£<\> 2 u - a£ t K + a£<f>tu - W^(f>xu ~ a£ 2 4> xu + (3£, 2 £ xx - a^t 2 + a£ t £ x ~ "CCxt = 0, 
3£,4> X u4 > xx + f3£4>4>xx - it4>xxx - 3^ xx (f> xx - £4> xx + a£4> x + 3£(/) xxu <t> x ^xxx^x 
- 2£,£x<f>x + it4>x + £,4>4>xxxu - £,4>4>xu + £</> xxx t - £<t> x t = o, 

3£,(f) u (t>xu - 6t,£x<f>xu - 3£ t (t>xu - ai 2 4>x - 3£, 2 4> xxu + - 3^ xx <f> u + 3£<f> xtu + K4>t 

+ Pttx<t> - PZt4> + ftixxx + Sttxtxx + 3^xx + 2t 2 tx - & X - 3^ xxt +6 = 0, 

S££ xx <t>xu + P£, 2 <Pxx + a£ t (f>x + 3££ x (f> xxu - 2ai4> u 4> x - Qi(t) 2 xu - a£(f> xt - ££ x £ xxx 

+ 2^4> xu + ££ xxx 4>u - 3(,4> u 4> xxu - 3^l x - 2f3i4>4> xu - 3£(f) xxtu + 2^ x 4> u - 2^ 2 X - £ 



(3.4i) 
(3.4ii) 
(3.4iii) 
(3.4iv) 
(3.4v) 

(3.4vi) 

(3.4vii) 



— a 



xxtu 
xu 



£t£x 



£,£.xx ~\~ S,^xxxt + 3e 



tYxxu 



0. 



(3.4viii) 
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These equations were calculated using the MACSYMA package symmgrp.max [15]. We then used the 
method of DGBs as outlined in [18,20] to solve this system. 

Case 3.1 a + (3 ^ 0. In this case it is straightforward to obtain from (3.4i-v), the condition 

£*^ 2 (a + /3)(3/?-2a) =0. 

The case 3/3 — 2a = leads to no solutions different from those obtained using the classical method. 
Subcase 3.1.1 £ x / 0. This is the generic case which has the solution 

£ = ( Kl x + K 2 )f(t), 4> = -«i/(t) (u - ^ + ^j^j + ^ 

where f(t) is an arbitrary differentiable function and K\ ^ and k 2 are arbitrary constants. These 
are equivalent to the infinitesimals (2.7) obtained using the classical method. 

Subcase 3.1.2 £ x = 0. In this case it is easy to obtain the condition 

^ 3 (/?-«)(« + /3) = 0. 

There are two subcases to consider. 

(i) a [3, £ x = 0. In this case the solution is 

£ = f(t), 4> = n 3 f(t) + 1/13, 

where f(t) is an arbitrary differentiable function and Ks is an arbitrary constant, which is equivalent 
to the infinitesimals (2.7) obtained using the classical method in the case when K\ = 0. 

(ii) a = f3, £ x = 0. In this case, we obtain the following DGB for £, cf) 

£ u = 0, i x = 0, 4> u = 0, 

£(f>xxxx ~ (£ + ^)<t>xx + P4>4>xx + P<f>l = 0, 

pe<Px - P&t + -6 = 0. 



Solving these yields 



where ( = x + /(i), /(i) is an arbitrary differentiable function and V(() satisfies 

V cc + f3V 2 -V = K4 ( + K 5 , (3.6) 

with k 4 and k 5 arbitrary constants. If re 4 / 0, then this equation is equivalent to the first Painleve 
equation [36] 

= 6y 2 + x, (3.7) 

otherwise it is equivalent to the Weierstrass elliptic function equation (2.14). 

Thus solving the characteristic equations (2.4) yields the nonclassical reduction 

u(x,t)=v(()+w(z)+t/f3, ( = x + f(t), z = x-f(t), (3.8) 

where f(t) is an arbitrary function and V(Q = and W(z) = w z satisfy 

V a + 13V 2 -V = -AC + fi u (3.9a) 
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and 

W zz + [3W 2 -W = Xz + fi 2 , (3.9b) 

respectively, where /ii and /i 2 are arbitrary constants and A is (effectively) a "separation" constant. 
If A / then these equations are equivalent to the first Painleve equation (3.7), whilst if A = 
then they are equivalent to the Weierstrass elliptic function equation (2.14). 

In particular, if A = (we set (3 = 1 without loss of generality) then equations (3.9) possess the 
special solutions 

V(C) = {6k 2 sech 2 (kiC) + 5 - 2k 2 } , W(z) = {6k 2 . sech 2 (« 2 z) + \ - 2k 2 .} , 

where k\ = | (1 + 4 / u :L ) 1 / 4 and k 2 = \ (1 + 4^ 2 ) 1 ^ 4 . Hence we obtain the exact solution of (1.8) 
given by 



u(x, t) = {6ki tanh {ki [x + f(t)]} + 6k 2 tanh {k 2 [x — f(t)]} 
+ x(l - 2k 2 - 2k 2 ) + 2/(t)(K 2 - k 2 ) + t}, 



(3.10) 



where f(t) is an arbitrary differentiable function. 

If Hi = jj.2 = then ki = k 2 = | (with /3 = 1) and (3.10) simplifies to 



u(x,t) = {3tanh{± [x + /(t)]} + 3tanh{± [x- f(t)]} + t} . (3.11) 

In Figure 1 we plot u x with u given by (3.11) for various choices of the arbitrary function f(t). This 
is one of the simplest, nontrivial family of solutions of (1.1) with a = (3{= 1), using this reduction. 
In Figure 1, f(t) is chosen so that f(t) ~ t + 1 , as t — > — oo, where to is a constant. Consequently 
all the solutions plotted in Figure 1 have a similar asymptotic behaviour as t — > — oo. However the 
asymptotic behaviours as t — > oo are radically different. 

In the special case when f(t) = ct, then choosing Ki = |(1 + 1/c) 1 / 2 and k 2 = |(1 — 1/c) 1 / 2 in 
(3.10) yields the two-soliton solution for (1.8) given by 



u(x, t) = —j= | y/c + T tanh ^(1 + 1/c) (x + ct) +y/c-l tanh ^(1 - 1/c) (x - ct) |. 

(3.12) 

This solution is of special interest since such two-soliton solutions are normally associated with so- 
called Lie-Backlund transformations (cf., [7]), whereas (3.12) has arisen from a Lie point symmetry, 
albeit nonclassical. Plots of (3.12) and its x-derivative are given in Figure 2(a) for c = 3. These 
should be compared with the plots in Figure 2(b) of the x-derivative of solution (4.6) below for the 
SWW-HS equation with ki = 2, k 2 = 1.7 and K\ = |, k 2 = | (so that A\ 2 = 1 in both cases). 

We stress that the "decoupling" of the nonclassical reduction (3.8) into a function of £ = x+f(t) 
and a function of z = x — f(t) occurs for the SWW equation (1.1) only in this special case when 
a = p. 

Case 3.2 a + (3 = 0. This case also leads to no solutions different from those obtained using the 
classical method in §2. 

To conclude this section we briefly consider the equations 

u xxxx + au x u xt + (3u t u xx - u xt - u xx = 0, (3.13) 
u xx tt + au x u xt + (3u t u xx - u xt - u xx = 0, (3.14) 

where a and (3 are arbitrary, nonzero, constants, which are variants of the SWW equation (1.1). 
Both these equations are thought to be nonintegrable, i.e. not solvable by inverse scattering, for all 
values of a and (3 since it is straightforward to show that neither satisfies the Painleve PDE test 
due to Weiss et al [55]. 
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Applying the classical Lie method to equations (3.13) and (3.14) yields the infinitesimals 



i = KlX + K 2 , T = (3«i + K 3 )t + K4, (/) 

£ = KiX + K2, T = — («l + «K 3 )t + K 4 , = (Ki + aK^U — K3X — 2(«i + OLKz)t/(3, 



(«1 + K 3 )ti 1 h K 5 , 

a p 



respectively, where k\, k 2 , ■ ■ ■ ,k$ are arbitrary constants. In contrast to the SWW equation (1.1), 
these symmetry groups are finite dimensional. 

If a = (3, then equations (3.13) and (3.14) both possess nonclassical reductions of the form 

u(x, t) = v(Q + w(z), Q = x + ct, z = x — ct 



where for (3.13), V(C) = an< ^ W(z) = w z satisfy 

V a - (1 + c)V + PcV 2 = -AC + Ai, 
W zz + (c - 1)W - pcW 2 =Xz + A 2 , 

and for (3.14), V(Q = and W(z) = w z satisfy 

cV cc - (1 + c)V + !5cV 2 = - AC + Ax, 
c 2 W zz + (c - 1)W - I3cW 2 = Xz + A 2 , 



(3.15i) 
(3.1Sii) 



(3.161) 
(3.16ii) 



with A, Ai and A 2 arbitrary constants. If A ^ 0, then equations (3.15,3.16) are solvable in terms of 
the first Painleve equation (3.7), whilst if A = then they are equivalent to the Weierstrass elliptic 
function equation (2.14). 

Setting A = Ai = A 2 = 0, and solving equations (3.15) and (3.16) yields the exact solutions of 
(3.13) and (3.14) with a = (3 given by 



u(x,t) = — {Vl + c tanh [|a/! + c(x + ct)] - y/l-c tanh [|a/I - c(x - ct)] } , (3.17) 



u(x, t) = ^ <J Vl + c tanh 



2c 



(x + ct) 



— Vl — c tanh 



2c 



(x — ct) 



(3.18) 



respectively, which are analogues of (3.12). Plots of the solutions (3.17) and (3.18) and their 
derivatives with respect to x are given in Figures 3 and 4, respectively. The solutions (3.17) and 
(3.18) are plotted in Figures 3(i) and 4(i) and look like the elastic interaction of a "kink" and an 
"anti-kink" solution. The x-derivatives of (3.17) and (3.18) are plotted in Figures 3(h) and 4(h) and 
look like the elastic interaction of two "soliton" solutions. These plots are very similar to the plots 
of (3.12) and its x-derivative given in Figure 2, however whereas (1.8) is completely integrable, both 
(3.13) and (3.14) are seemingly non-integrable. These solutions (3.17) and (3.18) are of particular 
interest since such solutions are normally associated with integrable equations, whereas they arise 
here for nonintegrable equations. Furthermore, like (3.12), these "two-soliton" solutions have also 
arisen from nonclassical reductions. As mentioned above, such solutions are normally associated 
with Lie-Backlund transformations (cf., [7]). 
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4 Soliton Solutions of the SWW-HS and SWW-AKNS Equations. 



4.1 The SWW-HS Equation. Exact solutions of the SWW-HS equation (1.8) can be obtained 
using the so-called singularity manifold method which uses truncated Painleve expansions [54,55]. 
If we seek a solution of (1.8) in the form 

6 Mx,t) u n 

HX,t) (3 <f>(x,t) 1 [ ' 

and then equate coefficients of powers of 4> to zero, we find that </>(x, t) satisfies the overdetermined 
system 

4>xxxt ~ <Pxx - 4>xt = 0, (4.2a) 

4>t4>xxx - ?>4>xt4>xx + 3<t>x4>xxt - 4>x(<t>x + 4>t) = o. (4.2b) 

(A DGB analysis of this system leads to some very complex expressions. Although it does yield 
some ODEs in x for <f> in the various subcases they appear difficult to solve.) Now suppose we seek 
a solution of these equations in the form 

4>(x, t) = ai exp {k\x + flit} + c*2 exp {k 2 x + fi 2 t} + ao, (4.3) 

where ao, «i, 0:2, fCi, K2, fii and fi2 are constants such that ao a i a 2 7^ 0. It is straightforward 
to show that equations (4.2) have a solution of the form (4.3) provided that fii = k\/(k\ — 1), 
fJ-2 = ^2/(^2 — 1) an< i K i an d ^2 satisfy the constraint 

k\ - K1K2 + kI = 3. (4.4) 

Thus we obtain the following exact solution of the SWW-HS equation (1.8) given by 

6 aiKiexp {kix + Kit/(Kl - I)} + a 2 K 2 exp \k 2 x + K 2 t/(Kl - 1)} 

u(x t) = — — (4.5) 

' /3 ao + «i exp {kix + Kit/(K^ — 1)} + a 2 exp {k 2 ^ + ^2^/(^2 — 1)} ' 

provided ki and K2 satisfy (4.4). It should be noted that (4.5) and (3.12) are fundamentally different 
solutions of the SWW-HS equation (1.8) as we shall now demonstrate. 
The general two-soliton solution of (1.8) is given by 

6_ Kiai exp (771) + k 2 os2 exp (772) + (ki + k 2 )Ai 2 exp (771 + t? 2 ) , 
(3 1 + Qi exp (771) + a 2 exp(77 2 ) + A i2 exp (771 +772) 

where 

77 1 = K lX + Kit/{Ki — 1), 7/2 = K 2 X + «2*/( K 2 ~ 1)' (4.6b) 

(«1 -K 2 ) 2 (Kf-KiK 2 + Ki-3) 

with Qi, CK2, «i and K2 arbitrary constants [34] (see also [43]). A straightforward method of 
obtaining this solution is to substitute (4.1) into SWW-HS equation (1.8) and then integrate once. 
This yields the bi-linear equation, 

4><PxXXt - <Pt4*XXX + 3(<f) xx (f) xt - (j> x (f> xxt ) - 4>{(f>XX + 4>xt) + 4>x((f>X + 4>t) = 0, (4.7) 

where we have set the function of integration to zero. Now, by seeking a solution of this equation, 
in the form 

(j)(x, t) = 1 + «i exp (771) + a 2 exp (772) + A12 exp (7/1 + 772) , 
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with r/i = K\x + fiit and r/ 2 = k 2 x + fi 2 t, it is easy to show that necessarily /ii = k,i/(k\ — 1), 
fj,2 = k 2 /(k 2 — 1) an d ^12 is given by (4.6c); this technique may be viewed as a variant of Hirota's 
bi-linear method [32]. Four plots of the x-derivative of (4.6), i.e. solutions of (1.6), are given in 
Figure 5 for (i), ki = 2, k 2 = 1.7, (ii), ki = §, k 2 = §, (hi), ki = 1.1, k 2 = (11 + 3\/93)/20 (so 
that A12 = 0), and (iv), «i = 0.8, k 2 = |(2 + 3\/7) (so that A 12 = 0). Plots 5(i) and (ii) illustrate 
"standard" two-soliton interaction whilst plots 5(iii) and (iv) are in the special case when Ai 2 = 0. 
Plots of the x-derivative of (4.6) for Figure 2(b) of the x-derivative of solution (4.6) below for the 
SWW-HS equation, i.e. solutions of (1.6), for k\ = 2, k 2 = 1.7 and ki = f , k 2 = § (so that A\ 2 = 1 
in both cases). 

The solution (4.5) is the special case of (4.6) with A 12 = 0; Hirota and Ito [33] refer to this as 
being the "resonant state" where either a single soliton splits into two solitons, see Figure 5 (hi), 
or two solitons fuse together after colliding with each other, see Figure 5(iv). On the other hand, 
(3.12) is the special case of (4.6) with Ai 2 = 1 where two solitons pass through each other with no 
phase shift as a consequence of the interaction. Thus whereas both solutions are asymptotically 
equivalent as t — > —00, they are qualitatively very different as t — > 00. This shows that the 
nonclassical method and the singular manifold method do not, in general, yield the same solution 
set. 

We remark that the SWW-HS equation (1.8) also possesses the solution 

, _ 6 Ki [«i exp (771) - a 2 exp (-r?i)] - K 2 £i2 sin (r? 2 + 5 ) g ^ 
(3 ati exp (771) + a 2 exp (—771) + B 12 cos (772 + 6 ) 



where 



(k| + k 2 2 - l)t 

Vl ~ KlX + (k? - 2 Kl + 4 + 1)(k? + 2«x +4 + iy [ } 

(4 + K 2 + 1)* / , , 

T72-K2X (At 2_ 2Ki+K 2 + 1)(K 2 +2Ki+K 2 + 1) ' i 4 -» C J 



b V 2 - ^ : - c;\ 1 :| V /2 ; (4.8d) 



^1/2 K i / 3 K i ^2 

K 2 V 3k 2 _ 4 + 3 

with cki, a 2 , Ki, k 2 and 5o arbitrary constants. Leble and Ustinov [38] show that solutions of the 
form (4.8) exist for several of PDEs that are solvable by inverse scattering through third order 
linear problems. Two plots of the x-derivative of (4.8), i.e. solutions of (1.6), are given in Figure 6 
for (i), K\ = 1, k 2 = 1.17 and (ii), K\ = 1, n 2 = 1.24. 

4.2 The SWW-AKNS Equation. The general two-soliton solution of the SWW-AKNS 
equation (1.7) is given by 

, t \_^_ exp (771) + K 2 a 2 exp (772) + (ki + k 2 )A 12 exp (771 + 772) ^ Q ^ 

(5 1 + ai exp (771) + a 2 exp (772) + ^12 exp (771 + 772) 

where 

771 = Kix + Kit/ (4 ~ !)> V2 = n 2 x + K 2 t/(4 - I), A 12 = aia 2 l I , (4.9b) 

\Ki + K 2 J 

with with ai, CK2 , Ki and k 2 arbitrary constants [34] (see also [43]). Two plots of the x-derivative of 
(4.9), i.e. solutions of (1.3), are given in Figure 7 for (i), ki = 2, k 2 = 1.7 and (ii), ki = |, k 2 = $. 

To apply the singularity manifold method to the SWW-AKNS equation (1.7), we seek a solution 
in the form 



Nonclassical symmetries and exact solutions of a shallow water wave equation 



13 



Equating coefficients of powers of (f> to zero, yields the overdetermined system 

fixxxxt — <Pxxx — 4>xxt = 0, (4.11a) 

<t>t<t>xxxx + Hx<f>xxxt - 1<t>xx<t>xxt - 3<t>x<t>xx - 4>t4>xx - 2<t>x4>xt = 0, (4.11b) 

ZMx^xxx - Hx<\>xx4>xt + 2<f>l<f>xxt - M4>1 + 4>lx) ~<f>l = 0- (4.11c) 



We note that this is a system of three equations in contrast to that for the SWW-HS equation 
(1.8) where <p(x,t) satisfies a system of two equations (4.2). It is easily shown that there exist no 
solutions of (4.11) in the form (4.3) if ao«ia2 / 0. Additionally substituting (4.10) into (1.7) yields 
a tri-linear equation for 4>(x,t), whereas the analogous equation for (1.8), i.e. (4.7), is bi-linear. 

5 Nonlinear Superposition for 2+1-dimensional Equations. 

In this section we show that nonlinear superposition of solutions can also occur for two 2 + 1- 
dimensional equations of the SWW-HS equation (1.8). Two 2 + 1-dimensional generalisations of 
the SWW equation (1.1) are the following equations 

Uyt + OLU x \l x y + (3llyU xx U XXX y 0, (5.1) 

U x t + &U x U X y + (3UyU xx U XXX y — 0. (^-2) 

We remark that these two equations differ only in the first term, which belongs to the linear part, 
yet their symmetries differ remarkably. Further both (5.1) and (5.2) reduce to the KdV equation 
(1.4) if y = x. 

Boiti, Leon, Manna and Pempinelli [13] developed an inverse scattering scheme to solve the 
Cauchy problem for (5.1) with a = (3, i.e. 

Uyt + U xxx y + (3Vj xx Uy + (3U x U x y 0, (5*^) 

which is a 2 + 1-dimensional generalisation of the SWW-HS equation (1.8), for initial data decaying 
sufficiently rapidly at infinity. This inverse scattering scheme is formulated as a nonlocal Riemann- 
Hilbert problem and involves a so-called "weak-Lax pair". We note that both (5.3) and SWW-HS 
(1.8) are reductions of the 3 + 1-dimensional equation 

Uyt + u XXX y + (5u xx Uy + [5u x u X y U xz 0, (5.4) 

which was introduced by Jimbo and Miwa [37] as the second equation in the so-called Kadomtsev- 
Petviashvili hierarchy, though (5.4) is not completely integrable in the usual sense (cf., [25]). 
Bogoyaviemskii [11,12] discusses the inverse scattering method of solution for (5.2) with a = 2(3, 
i.e., 

U X t + U XXX y + f3U XX Uy + 1(3U X U X y = 0, (^-^) 

which is a 2 + 1-dimensional generalisation of the SWW-AKNS equation (1.7); it should be noted 
that the 2 + 1-dimensional generalisation of (1.7) discussed by Gilson et aJ [28] is different to (5.2). 

It is routine to show that (5.1) satisfies the necessary conditions of the Painleve PDE test due 
to Weiss et aJ [55] to be completely integrable if and only if a = (3, i.e. when it reduces to (5.3), and 
(5.2) satisfies these necessary conditions if and only if a = 2/3, i.e. when it reduces to (5.5). This 
strongly suggests that (5.1) and (5.2) are solvable by inverse scattering only in these two special 
cases, both of which are known to be completely integrable ([13] and [11,12], respectively). 

5.1 Classical Symmetries. To apply the classical method to (5.1), we consider the one- 
parameter Lie group of infinitesimal transformations in (x, y, t, u) given by 

x = x + e£i(x,y,t,u) + 0(e 2 ), 
V = y + ^2(x,y,t,u) + 0(e 2 ), 
t = t + e£ 3 (x,y,t,u) + 0(e 2 ), 
u = u + e<j){x, y, t, u) + 0(e 2 ), 
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where e is the group parameter. Solving the resulting determining equations yields the infinitesimals 



x d^/i 

£l = KiX + /i(t), £ 2 =g(y), £ 3 = 3 Kl t + K 2 , (j) = -K 1 U J [ —+f2(t), 

a at 

if a ^ (5 and 



+ AW , 6 = 9 („), 6 -SAW, * = -f + ^ + ^f + / 3 W, 

if a = /3, where «i, k 2 and K3 are arbitrary constants and /i(t), fzif), fs{p) an d 9{y) are arbitrary 
differentiable functions [21] (classical symmetries in the case a = (5 are also discussed in [52]). 
Applying the classical Lie method to (5.2) yields the infinitesimals 

£l = KiX + /i(t), 6 = «2* + «3y + «4, £3 = (2Ki + K 3 )t + K 5 , 

, « d/l K2X 
= -KlU + -— + + / 2 (t), 

p di a 

if a 7^ 2/3 and 

6 = (^it + k 2 )x + /i(t), 6 = (2/«iy + k 3 )* + + K5, £3 = 2/^t 2 + (2re 2 + ^ 4 )t + k 6 , 
= -(k^ + k 2 )u + ^ ^KlX + ^\ + ^ + 

if a = 2/3, where aci, k 2 , . . . , kq are arbitrary constants and /i(t) and / 2 (t) are arbitrary differentiable 
functions [21]. 

5.2 Nonclassical Reductions. Applying the nonclassical method to (5.1) yields an additional 
reduction in the case when a = (3, when it reduces to (5.3), This nonclassical reduction is given by 

u(x,y,t)=v( V ,T)+w((,T)- r, = x+f(t)g(y), { = x-f(t)g(y), r = t, (5.6) 

where V(r), r) = v v and W((, r) = satisfy the variable-coefficient KdV equations 

V vriv - 2f3VV v — V T F{t)V + ^ |^ + 2F 2 (r)| + Q{t) = 0, (5.7) 

W ac - 2(3WW C -W T - F{t)W + ^ j ^ + 2F 2 (r) j + Q(r) = 0, (5.8) 

with F{t) = f(t)/f(t) and Q(t) an arbitrary function. It is easily shown that (5.7) and (5.8) satisfy 
the necessary conditions of the Painleve test due to Weiss et al [55] to be completely integrable and 
both equations can be transformed into the usual KdV equation (1.4). In particular, if in (5.7) and 
(5.8) we set f(t) = k, a constant, then we obtain solutions to (5.1) with a = (3 in terms of sums of 
solutions of the KdV equation (1.4), with arguments x ± g(y), where g(y ) is an arbitrary function; 
this is perhaps the simplest family of solutions to (5.3) found using this nonclassical reduction. 

Applying the nonclassical method to (5.2) yields an additional reduction in the case when a = j3, 
given by 

u{x,y,t) = v(t],t) +w((,t), 7] = x + ny, ( = x - ny, r = t, (5.9) 
where V(r), r) = v v and W((, r) = satisfy 

K%777j - 0kVV v -V t + Q(t) = 0, (5.10) 
/3kW ccc - f3KWW c -W T + Q{t) = 0, (5.11) 



with Q(t) an arbitrary function. 
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6 Discussion. 

In this paper we have discussed symmetry reductions and exact solutions for the shallow water 
wave equation (1.1). In particular, for the special case of (1.1) given by the SWW-HS equation 
(1.8), using the nonclassical symmetry reduction method originally proposed by Bluman and Cole 
[9], we obtained a family of solutions (3.10) which have a rich variety of qualitative behaviours. 
This is due to the freedom in the choice of the arbitrary function f(t). One can choose fi(t) and 
flit) such \f\{t) — f2{t)\ is exponentially small as t — ► — oo, yet fi(t) and fzif) are quite different 
as t — > — oo, so that as t — ► — oo the two solutions are essentially the same, yet as t — > oo they 
are radically different. In Figure 1 we show that by a judicious choice of f(t) we can exhibit a 
plethora of different solutions. We believe that these results suggest that solving the SWW-HS 
equation (1.8) numerically for initial conditions such as those in the solutions plotted in Figure 1 
could pose some fundamental difficulties. An exponentially small change in the initial data yields 
a fundamentally different solution as t — > oo. How can any numerical scheme in current use cope 
with such behaviour? 

The solution (3.10) appears to be a nonlinear superposition of solutions suggesting that the 
SWW-HS equation (1.8) may be linearisable through a transformation to a linear PDE, analogous 
to the linearisation of Burgers' equation 

ut=u xx +2uu x , (6.1) 

which is mapped to the linear heat equation through the Cole-Hopf transformation [22,35]. If so 
then the solution (3.10) could be viewed as an artefact of the fact that the SWW-HS equation 
(1.8) is linearisable. However the SWW-HS equation (1.8) can be expressed as the compatibility 
condition of the third order spectral problem (1.11,1.12). Further the associated scattering problem 
(1.11) is very similar to that for the Boussinesq equation which has been thoroughly studied by Deift 
et al [24]. This strongly suggests the SWW-HS equation (1.8) is solvable by inverse scattering. The 
spatial part of the inverse scattering formalism (1.11) only defines u up to an arbitrary additive 
function of t; this arbitrary function may be incorporated into u using the variable-coefficient 
transformation (2.15). Indeed the initial value problem for the SWW-HS equation (1.8) is not 
well-posed without the imposition of an additional constraint since if u(x,t) is a solution of (1.8) 
satisfying the inital condition u(x, 0) = <p(x), then so is 

u(x,t)=u{x,t) + [g(t)-t}/(3, x = x, t = g(t), (6.2) 

where g[t) is any differentiable function such that g(0) = 0. It appears likely that the inverse 
scattering formalism for the SWW-HS equation (1.8) will require that u(x,t) satisfies a constraint 
of the form 

/oo 
\u(x, t)\dx < oo, (6.3) 
-oo 

for all t. It is well-known that such constraints are required in the inverse scattering formalism for 
2 + 1-dimensional equations such as the Kadomtsev-Petviashvili equation (cf. , [6] , see also [1] and 
the references therein). 

Since the shallow water wave equation (1.1) is invariant under the variable-coefficient 
transformation (2.15) for all nonzero a and (3, one can take any solution of (1.1) and using (2.15) 
generate some interesting solutions. For example, one can apply the transformation (2.15) to the 
two-soliton solution (4.9) and thus generate some exotic solutions for the SWW-AKNS equation 
(1.7) analogous to the solution (3.10) of the SWW-HS equation (1.8). Consequently the above 
remarks about difficulties in solving (1.8) numerically apply to (1.1) in general. Additionally, the 
inverse scattering formalism for the SWW-AKNS equation (1.8) will probably require that u(x,t) 
satisfies a constraint such as (6.3). 

Recently Ablowitz, Schober and Herbst [5] have shown that the focusing nonlinear Schrodinger 
equation 

iu t + u xx + \u\ 2 u = 0, (6.4) 
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exhibits numerical chaos created by small errors on the order of roundoff. The results of Ablowitz 
et al [5] together with those given in this paper suggest that numerical analysts need to take care 
to ensure the accuracy of their programs. Numerical predictions of chaos may not always be what 
they seem! 
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Figure captions 

Figure 1. The solution (3.11) where 

(i) , /(*) = \t, 

(ii) , /(*) = |t-exp(it-l), 

(iii) , /(t) = |i+[l-tanht]sint, 

(iv) , it - exp (-t 2 /100) sint, 

(v) , /(*) = \t + 67T" 1 tan" 1 i, 

(vi) , f{t) = i^^T20, 

(vii) , f(t) = f(t) = \t(l - tanht) + |[2sin(|t) + 3](1 + tanht), 

(viii) , /(t) = ±t(l -tanht). 

Figure 2. (a) The solution (3.12) of the SWW-HS equation (1.8) and its ^-derivative for c = 3. 
Figure 2.(b) The x-derivative of solution (4.6) of the SWW-HS equation (1.8) for m = 2, k 2 = 1.7 
and Ki = | , K 2 = §• 

Figure 3. The solution (3.17) of equations (3.13) and its x-derivative for c=\. 

Figure 4. The solution (3.18) of equations (3.14) and its x-derivative for c = j. 

Figure 5. The solution (4.6) of the SWW-HS equation (1.8). 

(i) , K\ =2, k 2 = 1.7 

(ii) , Ki = §, K 2 = §, 

(iii) , ki = 1, k 2 = (11 + 3\/93 )/20 (A12 = 0), 

(iv) , ki =0.8, k 2 = 1(2 + 3^7) (A12 = 0). 

Figure 6. The solution (4.9) of the SWW-HS equation (1.8). 

Figure 7. Two soliton solutions (4.10) of the SWW-AKNS equation (1.7). 



